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Abstract 


Infrared dust bubbles play an important role in the study of star formation and the evolution of the interstellar 
medium. In this work, we study the infrared dust bubble N75 and the infrared dark cloud G38.93 mainly using the 
tracers C'8O0, HCO*, HNC and N>H* observed by the 30 m IRAM telescope. We also study the targets using data 
from large-scale surveys: GLIMPSE, MIPSGAL, GRS, NRAO VLA Sky Survey and Bolocam Galactic Plane 
Survey. We found that the C'8O emission is morphologically similar to the Spitzer IRAC 8.0 um emission. The 
1.1 mm cold dust emission of G38.93 shows an elongated structure from southwest to northeast. The ionized gas 
from G38.93 is surrounded by polycyclic aromatic hydrocarbon emission, which may be excited by radiation from 
G38.93. We found that the identified young stellar objects tend to cluster around G38.93 and are mostly in class II, 
with several class I cases distributed around N75, but no class II examples. We also found evidence of expanding 


feedback, which could have triggered star formation. 


Key words: ISM: bubbles — infrared: ISM — galaxies: star formation 


1. Introduction 


Churchwell et al. (2006) detected and cataloged about 600 
mid-infrared dust bubbles between longitudes —60° and +60°. 
Ionizing O- and/or B-type stars may play an important role in 
the formation of infrared (IR) dust bubbles (e.g., Zhang & 
Wang 2012, 2013; Zhang et al. 2013, 2016; Xu et al. 2014). 
These bubbles provide an important opportunity to study the 
interplay between the HII region and molecular clouds (e.g., 
Zhang & Wang 2012, 2013). Simpson et al. (2012) present a 
catalog of 5106 IR bubbles created by visual classification via 
the online citizen science website “The Milky Way Project.” 
Ultraviolet (UV) radiation from ionizing stars can heat dust and 
ionize gas, resulting in the formation of an expanding bubble 
shell (Watson et al. 2008). Massive star formation can be 
triggered by this process in the vicinity of shell clumps. 

Some bubbles are associated with infrared dark clouds 
(IRDCs). The IRDCs harbor various evolutionary stages, 
provide good targets to study the initial conditions of massive 
star formation prior to cloud collapse and are important to 
understand the early evolution of massive star-forming clumps 
(Rathborne et al. 2006; Beuther et al. 2007; Henning 
et al. 2010). Therefore, IRDCs play an important role in the 
study of star formation and the evolution of the interstellar 
medium (ISM) (e.g., Rathborne et al. 2005; Pillai et al. 2006; 
Beuther & Sridharan 2007; Zhang et al. 2017a). Additionally, 
bubbles associated with IRDCs will provide an important 


sample to study both triggered star formation and the early 
stage of high-mass star formation. The dense and cold nature of 
IRDCs may be more conducive to special chemical reactions 
than in other star formation regions (e.g., Zhang et al. 2017b, 
2019, 2020). For example, CO, a commonly used tracer of 
dynamics in cold and dense clouds (Hernandez & Tan 2011), 
can be affected by depletion of the gas phase abundance due to 
freezing on dust grains (Hernandez et al. 2011). Due to the 
freeze-resistance of some N- or D-containing molecules, N- 
and D-bearing molecules such as NH” are more suitable for 
tracing the dense and cold cores (e.g., Zhang et al. 2020). 

G38.93-0.355 (hereafter G38.93) is an IRDC associated with 
an IR dust bubble N75 (Churchwell et al. 2006). The distance 
of G38.93 was measured to be 2.7 kpc (Zhang et al. 2014). 
Bubble N75 and IRDC G38.93 are good targets to study the 
triggered star formation and high-mass star formation. The 
masers, including H2O masers (Forster & Caswell 1989, 1999), 
OH masers (Forster & Caswell 1989, 1999; Btaszkiewicz & 
Kus 2004; Szymczak & Gérard 2004; Caswell et al. 2013), and 
class II CHOH masers (Błaszkiewicz & Kus 2004; Szymczak 
& Gérard 2004), were detected in G38.93, suggesting massive 
star formation there. Xu et al. (2013) studied N75, N74 and 
their surrounding regions, however, they only used the 
optically thick CO as tracers to them. In this paper we use 
the new molecular spectra to perform a more detailed study of 
G38.93 to verify and complement the results of Xu 
et al. (2013). 
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In this paper we present a multiwavelength study in the 
bubble N75, mainly probing the ISM around N75 and 
searching for evidence of star formation. The observations 
and data reduction are described in Section 2; the results and 
data analysis are presented in Section 3; the discussion of the 
corresponding star formation is presented in Section 4. 
Section 5 is the conclusion. 


2. Observation and Data 


The adopted spectral data, covering the bubble N75 and the 
clump G38.93, are mainly from the IRAM 30m and FCRAO 
15m telescopes. They are used to study the dense and cold 
structure. In addition, the dust and continuum data are from 
Spitzer and the NRAO VLA Sky Survey (NVSS) archival data 
for studying the IR and HI background. Details about 
observations and data are presented in the following sections. 


2.1. IRAM 30m Observations 


The HCO* (1-0), HNC (1-0), NsH* (1-0) and C'8O data 
observations were made in 2018 December with the IRAM 
30m telescope? on Pico Veleta, Spain. The half-power 
beamwidth was between 29/0 and 23”5, the main beam 
efficiency (Beg) was between 81% and 78%, and the forward 
efficiency (Fett) of the IRAM 30 m telescope was between 95% 
and 94%. The relation between the main beam temperature 
(Twp) and the antenna temperature (Tx) is Typ = (Fore /Bere) X 
Tx. The on-the-fly mapping mode was used to scan each source 
in two orthogonal directions to reduce striping on the maps. 
The data reduction software is GILDAS.° 


2.2. CO Data 


The CO (J = 1 — 0) data are from the Boston University- 
Five College Radio Astronomy Observatory Galactic Ring 
Survey’ (GRS; Jackson et al. 2002). The data have a full width 
at half maximum (FWHM) of 46” and a velocity resolution of 
0.21 kms! at 110.2 GHz. The "CO data will be used to map 
our target and calculate the physical parameters. 


2.3. Dust and Continuum Data 


The IR data are from the Spitzer GLIMPSE survey 
(Benjamin et al. 2003) and the MIPSGAL survey (Rieke 
et al. 2004). The NVSS is a 1.4GHz continuum survey 
covering the entire sky north of —40° decl. NVSS was used to 
trace the ionizing gas (Condon et al. 1998) with a noise of 
about 0.45 mJy beam™'. We also used the 1.1 mm continuum 
data from the Bolocam Galactic Plane Survey (BGPS)® 
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(Aguirre et al. 2010) and the Bolocam 1.1 mm source catalog 
(Ginsburg et al. 2013) to trace the dense part of our target 
source. The effective FWHM of the 1.1 mm map is ~33”. 


3. Results and Analysis 
3.1. Dust Emission 


Figure 1 displays the '*CO emission (white), 1.1 mm 
emission (cyan) and 1.4GHz continuum emission (red) 
contours superimposed on the three-color image of 4.5 um 
(blue), 8.0 m (green) and 24 um (red). The Spitzer IRAC 8m 
emission may originate from polycyclic aromatic hydrocar- 
bons (PAHs), which are excited by the photodissociated 
region surrounding the HII region (Pomarés et al. 2009; 
Petriella et al. 2010). The 8 jm emission is extended at 
G38.93. The HII region N74, located to the southeast of N75, 
shows a ring-like shape. The 24 um emission, originating 
from hot dust, appears to be inside the bubble of N75 and 
N74. The 4.5 um emission indicates the location of the 
brightest stars. 

The 1.1 mm continuum emission is from cold dust. The cold 
dust emission of G38.93 shows an elongated structure from 
southwest to northeast. The structure of G38.93 is outlined by 
the 1.1 mm emission, which is comparable to the shape of the 
molecular emission contours in Figure 2. We used the NVSS 
1.4 GHz continuum emission overlaid on the three-color map 
to trace the ionized gas in the HII region. 


3.2. Integrated Intensity Maps 


We draw the integrated intensity maps of G38.93 and N75 in 
Figure 2. The BCO = 1-0) can be used to trace the structure 
of the clump associated with the HII region. We find that the 
composition of CO velocity, ranging from 35 to 42 km s™', is 
related to these H II regions. We present an integrated intensity 
map of '*CO(J= 1-0) in Figure 1, where two dense clumps, 
labeled with G38.93 and G38.95, can be clearly seen. 
Additionally, G38.93 is located at the northwest of N75, and 
has a compressed structure as seen from the CO distribution in 
Figure 1, suggesting a likely interaction between the N75 and 
G38.93. 

We also find that the HCO*, NH” and HNC emissions have 
similar morphologies to the 1.1 mm emission, while the co 
emission is more extended. The C'8O emission and the 8.0 um 
emission are very similar in morphology. In Figure 2 we see 
that the C'8O emission is weaker in the center of N75. The 
C'8O emission may be optically thick at such condition. 


3.3. Molecular Lines 


From the BGPS 1.1mm peak position of G38.93, we 
extracted individual co, HCO”, HNC and N>H* spectra. 
The average spectra of c'8o, HCOt, HNC and N>H* are 
displayed in Figure 3. These spectra show that the main 
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Figure 1. Three-color image toward H II region N75 and its surroundings with blue, green and red corresponding to 4.5 um (blue), GLIMPSE 8.0 jum (green; 
Benjamin et al. 2003) and MIPSGAL 24 jum (red; Carey et al. 2009), respectively. BGPS 1.1 mm emission is shown in blue contours, '*CO emission is in white 
contours and the NVSS 1.4 GHz radio continuum emission is in red contours. The white and blue contour levels are 10%-—100% of peak value by 10%. The symbols 
“E,” “a,” “x” and “x” mark the locations of CHOH class II maser, CH3OH class I masers, HO masers, and IRAS sources, respectively. 


location of the velocity components is in the range of 
33-45 kms '. We apply Gaussian fits to the spectra of C'%0, 
HCO*, HNC and NH”. The fitted results are summarized in 
Table 1. Then we derive a system velocity of ~38.98 km s~! 
using the C'5O line for clump G38.93. 


3.4. Velocity Distribution 


In Figure 4, we show the channel maps of c'80(1-0) to 
examine the detailed velocity distribution in the clump. We 
found that the system velocity of the bubble N75 is 
~40.5kms~'. Figure 5 presents the channel maps of the 
HCO” line. We can see that the 8.0 um grayscale shows the 
position and size of N74, N75 and G38.93. The HCO* emission 
extends over a velocity range of 34 < Visp < 42 km sl, while 
c'8O emission extends over a velocity range of 37.5< 
Visr < 41.5 km s™!. Based on the above velocity range, we 
make the position—velocity (PV) diagram of HCO* as displayed 
in Figure 8. Considering the PV diagram, we make the integral 
intensity maps of G38.93 as depicted in Figure 6. The blue 
contours indicate the blueshift velocity components and the red 
contours correspond to the redshift velocity components, 
suggesting the existence of outflow. 


Table 1 
Spectral Parameters of Molecular Cores at the G38.93 
Name Area Position Width Intensity 
(K km s7!) (km s~!) (km s~!) (K) 
c'o 2.729(0.056) 38.980(0.032)  3.031(0.069) 0.845 
HCO* 4.546(0.022) 40.659(0.007) 2.875(0.016) 1.485 
HNC 1.989(0.022) 40.610(0.012) 2.313(0.033) 0.807 
N>H* 1.851(0.017) 38.550(0.157) 3.850(0.157) 0.452 
8CO 13.884(0.381)  40.339(0.047) 3.47(0.108) 3.748 


3.5. Dust Temperature and Density 


In order to study the effects of the H II regions on G38.93, it is 
crucial to calculate the temperature and density. Herschel 
observations covering a wavelength range (70-500 um) can be 
used to construct the (Vy,) column density and dust temperature 
maps at 160, 250, 350 and 500 um. The reason we exclude 
70 um is that 70 um emission traces hotter components heated 
by UV emission. Using the Getsources algorithm (Men’ shchikov 
et al. 2010, 2012; Men’shchikov 2013), we created the column 
density and dust temperature maps of G38.93 by fitting the 
spectral energy distributions (SEDs) on a pixel-by-pixel basis 
(Palmeirim et al. 2013). 
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Figure 2. Integrated intensity maps of HCOt, C'80, HNC and N>H” lines on the Spitzer 8 um emission. The integrated velocity intervals are 33-45 km s~'. The 
HCO* contour levels start at 30 in steps of 30 (¢ = 0.1 Kkms_'). The beam size of each subfigure is indicated at the bottom-left contour. The straight line in the 


HCO* subfigure shows the cutting direction of the PV diagram in Figure 8. The unit of each color bar is in MJy sr_'. 


Figure 7 shows the Taust and Ny, maps. We found that the 
Taust distribution ranges from 10 to 30K, and the column 
density is about 10” cm? with a maximum column density of 
about 9.73 x 10° cm *. The temperature reaches minimum at 
the column density peak position of G38.93. Furthermore, the 
mean dust temperature is about 15.6 K, suggesting a relatively 
cold condition in this star formation region. 


3.6. Excitation Temperature, Column Density and Mass 
for G38.93 


Using the molecular lines, we calculate the physical 
parameters Tex, Na, m(H2) and Myre of G38.93 with a 
kinematic distance of 2.70 kpc (Roman-Duval et al. 2009). 
Assuming local thermodynamic equilibrium (LTE), we use the 
relatively optically thin '*CO(J = 1-0) emission to determine 


1 


the column densities and mass of the G38.93 clump by 
applying the following equation (Garden et al. 1991) 
(Tex + 0.89) 


N(@CO) = 4.57 x 10!83——* —_——_ 
exp(—5.29/Tex) 


f Tmpdv(cm~2), 
(1) 


where Tmb is the corrected main-beam temperature of G0, 
and Tex is the excitation temperature in K. Tex can be calculated 
by the following equation (Garden et al. 1991) 


5.29 


EF —. (2) 
In{1 + 5.29/(Tm (CO) + 0.89)] 


Tex 


The H, column density is calculated with the following 
equation (Frerking et al. 1982) 


N (H/N CO) = 5 x 10°. (3) 
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Figure 3. The spectra of C'8O, HNC, N»H* and HCO” lines at the clump G38.93. The green curves show the Gaussian fitting results. 
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Figure 4. Channel maps of the C'SO line over the velocity interval 35.5-42.5 kms’. The integrated velocity interval of each image is 1 km s~'. The central 


velocities are shown in each image. 


Then the mean volume densities of the Hə molecule are 
n(H2) = N (H2) /2R. (4) 


Based on the assumption of LTE and spherical geometry, the 
masses of the clouds are calculated as 


4 
Myre = pmyn(H2) x (tz), (5) 


where my is the mass of a hydrogen atom and u= 2.8 is the 
mean molecular weight considering the contributions of He and 
other heavy elements to the total mass. 

Finally, we derived the column density of the clump G38.93 
with M(H>) = 5.34 x 10°! cm~*. Other corresponding parameters 
are nH) = 1.16 x 10°cm™°, Rmin = 4/2, R=3.055pc, D= 
2.5 kpc, Tex = 6.953 K, N(3CO) = 1.068 x 10'° cm” and Mure = 
8 x 10* Mo. 
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Figure 5. The channel maps of the HCO* line over the velocity interval of 33-44 km s~". The integrated velocity interval of each slice is 1 km s 


velocities are given in each image. The blue contours range from 33 to 39 kms~' and the red contours range from 40 to 44 km s~". 


3.7. Outflow 


The PV diagram for G38.93 is visualized in Figure 8. The 
cutting direction passing through the selected dense core positions 
is indicated by the yellow line in the HCO” subfigure of Figure 2. 
The dense cores are the backbone of the G38.93 clump and allow 
us to study its dynamical state. The PV plot shows a clear velocity 
gradient relative to the system velocity of 39 km s_'. The broad 
and extended line wings of the HCO” in Figure 3 represent the 
blueshifted and redshifted velocity components respectively. The 
PV diagram in Figure 8 provides additional evidence for the 
outflow (Smith et al. 1997; Beuther et al. 2004). 

We selected the integrated area of the wings and determined 
the outflow intensities of the red and blue lobes according to 
the PV plots and channel maps. The contours of the integrated 
intensities of the HCO* line wing identified a molecular 
outflow (Figure 6). The dynamic timescale of the outflow is 
calculated using the equation tout = 9.78 x 10°R/ V(yr), where 
R in pc is the outflow size and Vin km s7" is the maximum 
flow velocity relative to the cloud system velocity. The mean 
dynamic timescale of the outflow is 5.1 x 10° yr. 


—! and central 


3.8. Expansion Age of the Bubble N75 


Using the 1.4 GHz radio continuum emission, the ionizing 
luminosity Nyy was computed via Mezger et al. (1974) 
with 


0.1 —0.45 2 
My = 4.76 x 108 4 ) (2) Sy) P|) s. 6 
GHz K Jy J\ kpc 


Here we adopted an effective electron temperature of 10° K, the 
distance 2.7 kpc and S$, = 23.7 mJy. We derived Ny = 2.14 x 
10f6s7!. Assuming the HI region is expanding into a 
homogeneous medium, the dynamical age was estimated by 
the method of Dyson & Williams (1980) 


4/3 —0.25 —0.5 
Ryu My ( nj ) 
tuu = 7.2 x 104] — — T. 
ii ( pe (a s7! 103 cm~? j 


(7) 


The age of the H N region G38.93 is 1.38 x 10° yr. 
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Figure 6. Integrated intensity map of the HCO* outflow (red and blue contours) superimposed on the Spitzer 8 jm emission. The red and blue contours are 30%, .. 
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100% of the peak. The red squares and triangle represent the IRAS sources and the Class II methanol maser (Deharveng et al. 2010) respectively. The beam size is 
given in the lower left corner. The symbols “9,” “A,” “x,” “W,” “a” and “x” indicate the location of an OH maser, IRAS, CH3OH masers, HO masers and the peak 


positions, respectively. The green cross marks a BGPS source. 
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Figure 7. Left panel: the BGPS 1.1 mm continuum emission (red contours) overlaid on the Herschel temperature map. Right panel: the BGPS 1.1 mm continuum 
emission (red contours) superimposed on the Herschel column density. Both maps were constructed using the SED pixel-by-pixel fit. 


4. Discussion 
4.1. Star Formation Activity 


The molecular outflows (see Figure 6) and masers such 
as H20, OH, and Class II CH30H (Caswell et al. 1995) 
were detected surrounding the bubble N75 and G38.93. This 
indicates that such a star formation system is an active 
star formation region. Furthermore, we also found many 
young stellar objects (YSOs) associated with this star 


formation system, so that we could derive their spatial 
distribution, ages, and masses. Based on the criteria of Allen 
et al. (2004), we selected near-infrared sources with 3.6, 4.5, 
5.8, and 8.0 um emission in our survey region from the 
GLIMPSE point source catalog. Figure 9 shows the [5.8]- 
[8.0] versus [3.6]-[4.5] color-color diagram. Only sources 
detectable in four Spitzer IRAC bands were considered. In 
Figure 9, class I YSOs are protostars with circumstellar 
envelopes, class II YSOs are disk-dominated objects and 
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Figure 8. PV diagram of HCO* along the straight line of the clump G38.93 in 
Figure 2. The contour levels start from 5ø to the peak integrated intensity of 50 
(o = 0.146 K). The gray bar represents the intensity in Kelvin. 


class II cases are main sequence or giant stars. From Figure 10, 
YSOs tend to cluster around G38.93 and mostly in class II, 
with several class I instances distributed around N75, but 
no class II cases. The ages of class I YSOs are ~10° yr, 
while the lifetimes of class II YSOs are ~10° yr (Andre & 
Montmerle 1994). 
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Figure 9. GLIMPSE color-color diagram [5.8]-[8.0] vs. [3.6]-[4.5] for 
sources in and around G38.93. Class I, II and III sources are indicated with 
different color areas. The classification as Class I, II, or IN suggests the stellar 
evolutionary stage as defined by Allen et al. (2004). 


4.2. Triggered Star Formation 


Figure 10 shows an RGB image with 4.5 wm, 8 um and 
24 um in blue, green, and red, respectively. The positions of 
the candidate YSOs are shown in Figure 10, where Class I and 
Class II YSOs are indicated by black and red dots, 


032450 


033650 


Figure 10. Three-color image with IRAC 4.5 jum, 8 wm and MIPS 24 um in blue, green, and red respectively. Class I and II YSOs are marked with black and red dots 
respectively. White outlines are '*CO emission. The square indicates the area for the identified YSOs. 
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respectively. The YSOs are mainly located in the compressed 
regions between bubble N75 and G38.93, suggesting that 
expanding feedback could have triggered star formation. 
Four Class I stars are located on the compact shell of bubble 
N75, which may indicate triggered star formation. The age of 
N75 is larger than that of those Class I stars, suggesting 
that there is enough time for this trigger process. Class II stars 
are concentrated at the east of G38.93 with strong 8 um 
emission. The 8 um emission distribution has a cocoon-like 
structure, which could be a compact HII region ionized by 
massive star formation. It is probably that the expansion of 
such a compact H II region has triggered the formation of class 
II stars, the masers and the outflows. We note that only a few 
YSOs are consistent with the bright bubble edges. It is 
possible that some YSOs are simply hidden in the high IR 
background, and the YSOs are mainly located in the 
compressed regions between the bubble N75 and G38.93. 
This may suggest that expanding feedback could have 
triggered star formation. 


5. Conclusions 


In this work, we studied bubble N75 with c!8o, HCO’, 
HNC, and NoH", and performed a multiwavelength study 
around G38.93. We found morphologically the C'8O emission 
correlates well with the Spitzer IRAC 8.0 wm emission. The 
east of G38.93 shows an arc-like shape, with a small area of 
higher PAH emission in the center. The 1.1 mm emission, 
tracing the cold dust, outlines the shape of G38.93, which is 
similar to the shape of the molecular emission. The cold dust 
emission of G38.93 exhibits an elongated structure from 
southwest to northeast. The 1.4GHz continuum emission 
shows the ionized gas in G38.93. The ionized gas from G39.93 
is surrounded by the PAH emission, which may be excited by 
radiation from G38.93. We found that the identified YSOs tend 
to cluster around G38.93 and mostly in class II, with several 
class I cases distributed around N75, but no class II examples. 
We also found evidence of expanding feedback, which could 
have triggered star formation. 
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